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I. INTRODUCTION 

Magnetic multilayers and spin valves exhibit a rich and 
extremely useful array of magneto-electronic phenomena, 
including magnetoresistance, interlayer exchange cou- 
pling, and spin transfer torque. Spin transfer torque is 
an effect in which the application of a current can change 
the relative orientation of the layers' magnetization. The 
presence of spin transfer torques are understood as a 
result of noncquilibrium electron current flow, or more 
precisely, a spatially varying spin current. Spin currents 
describe the flow of electron spin, and are tensor objects 
with two indices: one labeling the direction of flow in real 
space, the other labeling the vector component of spin. 
Interlayer exchange coupling is an indirect exchange in- 
teraction between magnetic layers, like the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction, with equilib- 
rium electrons in the nonmagnetic spacer playing the me- 
diating role. Interlayer exchange coupling is closely re- 
lated to spin transfer torques, as we discuss below. In 
this paper, we elucidate the relationship between inter- 
layer exchange coupling and spin transfer torques, specif- 
ically how the interlayer exchange coupling depends on 
the applied bias. 

Interlayer exchange coupling gives rise to an ener- 
getic preference for parallel or anti-parallel alignment of 
the layers' magnetization, and is given approximately as 
E ~ — JM 1 ■ M2 , where Mi are the directions of the mag- 
netizations. It shows RKKY-like characteristics, such 
as an oscillatory dependence on spacer layer thickness, 
with a period determined by the spacer material Fermi 
wave vector. This behavior is well establishedii^ by the 
close comparison between the experimental^!^ and theo- 
retical results^iS'i. The interlayer exchange coupling can 
be computed using the general formalism for calculating 
magnetic exchange interactions in bulk ferromagnets^. 
Slonczewski has pointed out that interlayer exchange 
coupling can also be understood as a consequence of con- 
servation of total spin angular momentum and the pres- 
ence of spatially varying spin currents in equilibrium (by 



equilibrium we mean no net charge current flow)^. In a 
spin valve, these equilibrium spin currents' spin vector 
is perpendicular to the magnetization of both magnetic 
layers, or out of the plane spanned by the two magne- 
tizations. This out-of-plane spin current results in an 
out-of-plane torque, denoted by r_L, which is responsible 
for the interlayer exchange coupling. Its magnitude is 
related to the interlayer exchange coupling, and is given 
by = J\Mi X M2I. 

From the perspective of magnetic exchange interac- 
tions as arising from the presence of spatially varying 
spin currents, it is a small conceptual step to deduce the 
presence of spin transfer torque when there is a nonzero 
net charge current. Indeed, typical spin transfer torques 
are the result of a net spin current with spin vector in 
the plane spanned by the two magnetizations^ i^°'^^ , lead- 
ing to torques in the plane, denoted by Fy. Frequently, 
the equilibrium torque is referred to as the interlayer ex- 
change coupling and the additional torque due to cur- 
rent is referred to as the spin-transfer torque. We find 
it conceptually useful to refer to the in-plane torque as 
the spin transfer torque and the equilibrium plus the 
bias-dependent out-of-plane torque as the interlayer ex- 
change coupling. (Note others have referred to the bias- 
dependent part of r_L as the "field-like spin transfer 
term" , as its form is the same as the torque from an 
applied field.) The physical difference between interlayer 
exchange coupling and spin transfer torque is then due 
simply to the different spin vector components of the spin 
current (see Fig. [T|). We use the terms interlayer ex- 
change coupling and out-of-plane torque Fj^ interchange- 
ably. 

The fact that Fj^ and F|| originate from different com- 
ponents of the spin current's spin vector imply further 
qualitative differences between the torques. The most ob- 
vious difference is that Fj^ is present in equilibrium, while 
F|| is only present in noncquilibrium systems. This differ- 
ence can be understood on general symmetry grounds^^. 
Another important difference is the relative magnitude 
of F^ and Fy. For current densities attained in mag- 
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netization switching experiments (10 A/m ), Fy ^ T± 
for metallic systems, which can be understood on gen- 
eral grounds and is discussed further in Section (jlll Ap . 
Finally, there is generally a qualitative difference in the 
dependence of the in-plane and out-of-plane torques on 
the applied bias. (We note that this different bias- 
dependence of Fj_ and Fy renders comparisons of their 
magnitude only partially meaningful.) This difference 
in bias-dependence is the motivation for this work, as 
there are inconsistent experimental results on the out- 
of-plane torque bias-dependence. For tunneling systems, 
a dependence of the interlayer exchange coupling on the 
measured bias potential has been reported as linearis, or 
quadratioi^, while for metallic systems, experiments have 
shown that interlayer exchange coupling is nearly inde- 
pendent of bias^"*, while others are interpreted as indicat- 
ing a linear dependence with substantial slope^^. Theo- 
retical work has implied a linear relation of the interlayer 
exchange coupling with bias for metallic systems^^iH, or 
posited a quadratic dependenc o^^'^^ . Ref. [2d| proves on 
symmetry grounds a quadratic dependence for symmetric 
systems with semi-infinite ferromagnetic leads. 
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FIG. 1: Spin valve showing the different components of 
torques on a magnetic layer. Bold arrows indicate the direc- 
tion of magnetization in the two ferromagnetic layers (light 
grey) surrounded by nonmagnetic material (dark gray). The 
assumed potential drop is shown above the multilayer: the 
dark curve shows the linear potential drop used for the nu- 
merical results, while the dashed step-like function is the po- 
tential drop assumed in analytic formulas. Contributions to 
the torque are found from states emanating from L and R 
lead separately. 



In light of these disparate results, this work consid- 
ers more carefully the bias dependence of interlayer ex- 
change coupling (or Fj^). The question of how magnetic 
exchange interactions are altered in the presence of cur- 
rent flow has been considered before in Ref. [2l|, which 
studied how the RKKY interaction between two spins 
changes due to a current-carrying electron distribution, 
and in Ref. [H, which finds a bias-dependent exchange 
interaction proportional to the ferromagnetic layer thick- 
ness. Heidc proposed that effects now understood to be 
the result of spin transfer torques were a manifestation 
of current-altered exchange interactions22i. We add to 
these works in the hope of addressing recent experimen- 
tal questions, and to express bias-dependent exchange 
interactions in a language familiar to previous studies 
of equilibrium exchange interactions and interlayer ex- 
change coupling. Our general results are valid for both 
multilayers and bulk ferromagnets in the regime of bal- 
listic transport. 

We restrict our attention here to spin valves. We find 
that the bias dependence of interlayer exchange coupling 
can exhibit a wide array of behaviors, so we focus here on 
some special cases in which the physics is clear and which 
are relevant for realistic systems. We find for metallic sys- 
tems in the limit of weak magnetic potentials and large 
spacer thickness that the interlayer exchange coupling 
shows oscillatory dependence on bias; however for exper- 
imentally attainable bias potentials, the change in inter- 
layer exchange coupling is negligible. In the symmetric 
half-metallic metal and tunneling limits, the interlayer 
exchange coupling depends quadratically on the bias. We 
derive approximate expressions for the prefactor of this 
quadratic term and compare them to exact numerical re- 
sults. 



II. METHOD 

The approach we adopt is relevant for systems in the 
ballistic limit. In the spirit of Landauer's approach to 
transport in mesoscopic systems, we suppose the system 
is connected to two reservoirs (to the left (L) and right 
(R)) in separate thermal equilibria, with chemical poten- 
tials /iL and /iR . A bias voltage eV across the system 
is represented by a difference eV between the reservoir 
chemical potentials. States in the system can be classi- 
fied as emanating either from the L or R lead (our formal- 
ism omits contributions from bound states with energies 
within the bias window). We find expressions for the in- 
terlayer exchange coupling due to states emanating from 
the L and R lead separately. This permits us to find the 
interlayer exchange coupling when L and R leads have 
different chemical potentials (that is, when there is an 
applied bias). 

We utilize a noncquilibrium Green's func- 
tion approach, in which the presence of the 
Left/Right lead is accounted for with a retarded 

(advanced) self-energy Tr^j^^{E), and observ- 
ables are expressed in terms of Green's functions 

G'-^'^HE) = (^E - H -1:1^''\e) -i:'^''\E)y\ Our 

approach generalizes the expressions of Liechtenstein 
et al.^ to find magnetic exchange interactions for 
noncquilibrium systems. 

When the spins of some subset of orbitals {i} are ro- 
tated by an angle 6 from their coUinear ground state 
orientation, there will be exchange torques on these or- 
bitals of the form J = Ji sin(6') -I- J3sin^(6') -I- ... . As- 
suming /iL > ^J■R, we state here the final result for ex- 
change torque coefficient Ji (we defer the derivation, and 
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more details regarding nonequilibrium Green's functions, 
to appendix A): 



1 (-MR 

Jl = -Im / Tr [A {G\ - G\) - AGf Aq] dE 

^ J ~OQ 



Re Tr 

J fJ^R. 



a(g< 



(AGf AG<^ + AG<^AGJ 



dE 
(1) 



where A = iJ | — is the spin-dependent Hamiltonian, 
Gf — T7-G!l (Ef — S?) G", and the traces are over or- 
bitals {i} that have rotated spins. For V = 0, the limits 
on the 2nd integral are identical so that its contribution 
vanishes. In this case Eq. IH) reduces to the well estab- 
lished expression for equilibrium exchange interaction of 
Ref. [^. For V ^ 0, in addition to contributions from the 
second term of Eq. ([T]) , the bias potential will also alter 
the Greens functions. Eq. ([T]) is applicable for both bulk 
ferromagnets and magnetic multilayers. (In our previous 
ab initio work on spin transfer torques in Co-Cu spin 
valves, we presented only the contribution from the sec- 
ond term of Eq. ^^^.) 

As an application of Eq. ([1]), we consider a one- 
dimensional, single band tight-binding model of two fi- 
nite ferromagnetic layers with riFi,nF2 sites, separated 
by a nonmagnetic spacer with ns sites, with nonmagnetic 
leads. The model is parameterized by the spin-splitting 
of the ferromagnetic (FM) layers (assumed equal) A, a 
scalar potential in the spacer layer [/, the Fermi energy 
Ey, and the number of layer sites npi, nF2, Jt-s- n-s de- 
termines the spacer layer width D: D = a{ns + 1) for a 
lattice spacing a. Energies are presented in dimension- 
less form, scaled by the tight-binding hopping parameter 
t. We focus here on some limiting cases in which the 
physics is clear, and which are relevant for realistic sys- 
tems. In all of our numerical calculations we assume 
that the bias potential drop is linear and occurs over the 
spacer layer (see Fig. ([!])), although in our hmiting ana- 
lytic results we neglect this linear potential and assume 
a step-like potential. 



III. RESULTS 

A. Small confinement, large spacer 

In the spirit of earlier work on interlayer exchange cou- 
pling, we first consider weak spin-dependent potential in 
the ferromagnet V^^^, and large spacer thickness D. In 
this limit, the one-dimensional, zero bias interlayer ex- 
change coupling takes the well-established form^i^ii: 



J 



where '''1(2) reflection off of the spin- up (down) 

potential of layer 1 (2), and k{EF) is the Fermi wave- 
vector of the spacer layer, vp is the Fermi velocity, given 
by vp = dE/hdk. (Note that 'f'^i(^2) includes the thickness 
dependence of the ferromagnetic layers.) The interpreta- 
tion of Eq. ^ is that to lowest order in r, the two layers 
are magnetically coupled together via itinerant states in 
the spacer through spin-dependent reflection of spacer 
states on both layer 1 and 2, with the mediating state 
picking up a phase factor 2k{Ep)D for the round trip 
between the layers. Eq. ([2]) has been successful in de- 
scribing the spacer material and thickness dependence of 
interlayer exchange coupling in experimenta^i^i^. 

At zero bias, the interlayer exchange coupling on layers 
1 and 2 is necessarily equal and opposite, because there 
is no external source of angular momentum. Therefore 
the choice of rotating layer 1 or 2 to flnd the interlayer 
exchange coupling is arbitrary. Under bias, the interlayer 
exchange coupling on layers 1 and 2 can be different be- 
cause the reservoirs supply angular momentum when out 
of equilibrium. This implies that the interlayer exchange 
torque is not the same on layers 1 and 2. Here we consider 
the interlayer exchange torque on layer 1. As described 
in the previous section, we flnd the torque Jl/r due to 
electrons from the L/R lead separately. The calculation 

of Jl/r 
result: 



can be found in Appendix B, here we give the 
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We flnd that that all of the out-of-plane torque on layer 
1 is due to states emanating from L, while states from R 
contribute nothing. This can be understood in terms of 
paths of states from L or R considered to second order in 
the reflectivity in r. Fig. ([2]) shows the relevant path for 
states from the L and R leads. A state from L transmits 
through 1, reflects off 2 (so that it now "carries" informa- 
tion about 2), comes back and reflects off 1, effectively 
"communicating" that information to 1, thereby coupling 
the two layers. For states from R, the two-reflection path 
transmits through 2 (so is initially oblivious of it), reflects 
off 1, comes back and reflects off 2. A final reflection 
off 1 would communicate spin information from 2, but 
this would be a third order process. There is therefore 
no second order contribution to the interlayer exchange 
coupling on layer 1 with this path. There are, of course, 
contributions to both Jl and Jr that are higher order in 
the reflectivities. However, in the limit that the reflectiv- 
ities are much less than one, the additional contributions 
are negligible. 

We assume that the dependence of J on Ep is 
dominated by the exponential factor. (That is, we 
assume D (dk/dE) is much larger than dr/dE and 
{h/ JD){dvp/dE), which is valid in the large D limit). 
The bias dependence can then easily be found: 

hvp 
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FIG. 2: Paths for states emanating from L and R leads 
for interlayer exchange couphng to lowest order in the spin- 
dependent reflectivity Ar. 



where eVc is the bias potential difference between left 
lead and spacer layer, and the wave vector fc is a func- 
tion of the argument Ep + gVl- The interlayer exchange 
coupling is an oscillatory function of the bias V, with 
period 2'Khvp/D. Fig ^ shows numerical evaluation of 
Eq. ^ for the 1-dimensional tight binding model, and 
we find excellent agreement with the approximate J{V) 
of Eq. (|4]) in the large D regime. We note that the phase 
of the oscillation is unconstrained, and is set by model 
parameters, most experimentally relevantly by D. 



where e is the electron charge, and |Jo| is the max- 
imum equilibrium value (in the 1-d model, |Jol = 
?iz;F|AriAr2|/27rD), which in this approximation is also 
related to the value of the coupling at the extremal wave 
vector, and includes prefactors such as the curvature of 
the Fermi surface at the extremal point, (see Refs. 00 
for example) , and is the Fermi velocity of the extremal 
wave vector. In the case of multiple extremal wave vec- 
tors, Eq. ([5]) should be a sum over all of them. For a Cu 
spacer with D = IQ nm, the change in J is | Jo|0.971 V. 
Assuming the maximum voltage on a metallic multilayer 
to be on the order of 10"''' V, this leads to a maximum 
change on the order of 10~^| Jo| (typical values for Jq are 
0.24 mJ/m^f21). We therefore conclude that in metallic 
systems, the bias-induced change in interlayer exchange 
coupling is negligible and not measurable. 



The oscillatory dependence of the exchange coupling 
(or out-of-plane torque) on bias may be surprising in light 
of the general argument for quadratic dependence in sym- 
metric junctions given in Ref. 20. It is important to note 
that this argument assumes semi-infinite ferromagnetic 
layers, and therefore does not apply to the system under 
consideration here. 
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FIG. 3: {J/t ■ D/a) vs. bias potential for several spacer layer 
thicknesses, npi = np2 = 7, A = O.lt, (7 = 0, Ef ^ -0.78f 
for our 1-d model. Legend shows different spacer widths in 
units of the lattice spacing a. 



The success of Eq. ([2]) in describing experimental data 
encourages us to believe that Eq. ^ is also experimen- 
tally relevant for metallic systems. Making a stationary 
phase approximation and considering only extremal wave 
vectors of the Fermi surface, the maximal possible slope 
of J{V) near bias is: 



It is instructive at this point to contrast the physics 
of current-altered interlayer exchange coupling and spin 
transfer torques. We emphasize that the interlayer ex- 
change coupling properties are dominated by spacer state 
properties (such as D, fcp, and the curvature of the 
Fermi surface around extremal points); in this limit, a 
bias changes interlayer exchange coupling to the extent 
it changes these spacer state properties. The sensitivity 
of interlayer exchange coupling to spacer state properties 
can be motivated by considering the interlayer exchange 
coupling as resulting from out-of-plane spin currents. For 
layers 1 and 2 in the x and z directions respectively, there 
is no obvious natural preference for the out-of-plane spin 
current to be in the positive or negative y direction. In- 
deed, the values of out-of-plane spin current for different 
states are distributed around 0, and mostly cancel. Only 
around extremal points is there a net contribution (this 
is why expressions for interlayer exchange coupling in- 
volve quantities evaluated only at E-p and at extremal 
wave vectors), and the final value is in some sense acci- 
dental. On the other hand, the in-plane torque acting on 
layer 2 has a strong preference to point in the x direc- 
tion. As for the out-of-plane torque, the specific values of 
in-plane torque from different states will vary, but they 
are distributed now about a nonzero value. For typical 
metallic systems such as Co-Cu at a current density of 
J — lO^A/cm^, the magnitude of the out-of-plane torque 
is 5 % that of the in-plane, reflecting this qualitative dif- 
ference in the spin current vector. This point has been 
made by Bauer et al. in the language of the mixing con- 
ductance G^-'-. There the effect is manifest in a small 
imaginary part of G^-^^^. 
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B. Half- metal 



In this section we take the hmit of a metaUic junction 
with half-metalhc FM layers. We choose this system be- 
cause it most easily demonstrates the opposite limit of 
the small confinement case. We find that in the half- 
metallic case, electrons from the L and R leads contribute 
equally to the interlayer exchange coupling on each layer 
(whereas in the small confinement case, electrons from 
only one lead contributed to the interlayer exchange cou- 
pling). Here we assume the potential for majority chan- 
nel is flat throughout the magnetic and nonmagnetic lay- 
ers, and the minority potential is infinite in the ferromag- 
net. The out-of-plane torque coefficient (near antiparallel 
alignment, or 6* « tt) then takes the simple form: 



sin [2k (E) D] dE , 



(6) 



This expression can be found by computing the out-of- 
plane spin current in the spacer layer. Integrating over 
all energies leads to the exact result for the interlayer 
exchange coupling: 



Ji.,K — 



2vYh{D/a^) cos{2kFD) + Ep sm{2kFD) 



A{D/a) 



1 



(7) 



The bias dependence is given by simply shifting the 
energies of L and R electrons by +V/2 and —V/2, re- 
spectively, leading to: 



JiV) - Jl LBf + - 



J(0) 



Jr 

Id'^JiEp) (V 



dE^ 



(8) 



Fig. [5] shows numeric results of the bias dependent in- 
terlayer exchange coupling for a number of spacer widths. 
We find these results correspond well with Eq. ([5]) . The 
curves for interlayer exchange coupling are again peri- 
odic in V, with the same period 2'Khvp/D as the previous 
small r case. In contrast to the small r case, the phase 
of the oscillations are constrained so that all curves are 
cos-like, which leads to a quadratic dependence of the 
exchange torque at small V. 

The half-metallic system is emblematic of the general 
fact that if states from the L and R leads contribute 
equally to the interlayer exchange coupling, and if there 
is a symmetric bias potential drop, then the bias de- 
pendence of the out-of-plane torque must be quadratic. 
Mathematically this quadratic dependence is seen imme- 
diately from Eq. ([5]) . A physical description is that when 
a bias eV is applied, we gain contributions from L states 
at energies above the original Fermi energy (from Ep to 
Ef + ^ ), while we lose contributions from a lack of R 
states at energies below the original Fermi energy (from 
i?F — 3" to Ep) (see Fig. dS])). This gain and loss cancel 
each other to linear order in eV, and only a quadratic 
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FIG. 4: {J/t ■ D/a) vs. bias potential for several spacer layer 
thicknesses, npi = nF2 = 7, -Ef = — 0.78f. FM are half 
metals and results are for 6 near tt. Legend shows different 
spacer widths in units of the lattice spacing a. 
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FIG. 5: Schematic of how bias potential changes the contribu- 
tions from the L and R leads to the total interlayer exchange 
coupling. At zero bias both L and R are contribute equally 
up to E-p, while a symmetric bias drop of eV opens a window 
from E-p — eV/2 to E-p -I- eV/2 with only one lead contributing. 
This leads to a quadratic dependence of the coupling on the 
bias. 



dependence remains^^. It should be noted that in mak- 
ing a constant shift in energies of the left and right leads 
relative to the spacer layer, there are also changes in the 
contribution from the bottom of the bands. However near 
the bottom of the band, the contributions to r_L oscillate 
rapidly and don't contribute (due to a large so that 
shifts in this region are unimportant. The quadratic de- 
pendence of the out-of-plane torque is also discussed in 
Refs. IT9H20I for the case of semi-infinite ferromagnetic 
layers. 
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C. Tunneling 

We next consider the case of a tunneling barrier with 
half-metaUic ferromagnetic layers. Such a model system 
is relevant for the states at the Brouillon zone center of 
Fe-MgO-Fe. These states make the dominant contribu- 
tion to the transport^l, spin torques^Si, and interlayer 
exchange coupling^^. An exact expression for the out-of- 
plane torque is given in Ref. [ol for a trilayer (FM leads 
separated by a spacer); for our 5-layer geometry, F^ is 
generally a complicated function of model parameters. 
We find that in the symmetric system (npi = 71^2), elec- 
trons from both L and R leads contribute equally to the 
interlayer exchange coupling, so that there is a quadratic 
dependence of the interlayer exchange coupling on bias 
voltage, as discussed in the previous section. We make a 
simplifying ansatz here for Jl/r and compare its predic- 
tion to the prefactor with the exact numerical results. 
We find that our approximate, analytic form shows qual- 
itative agreement with the numerics. 

We suppose that the energy dependence of the out-of- 
plane torque is dominated by the tunneling exponential: 



(9) 



where k is the decay constant (or imaginary wave- 
vector) in the spacer. In our model, Ka = 
Im [cos^^ {{Ef — U) /2t)] . ^ is a complicated function 
of model parameters, and we have found that compared 
to the exponential factor, A depends weakly on energy; 
we therefore omit its energy dependence going forward. 
To find the interlayer exchange coupling, wc assume ener- 
gies near the Fermi level make the dominant contribution 
to the total, leading to: 



Atsinh(KFa) _2kfD 
Jl,r = jj, e 



(10) 



where D' = D/a. The bias-dependent interlayer ex- 
change coupling can be written in the form: 



JiV) = J(0) (l + BV^) . 
A straightforward calculation leads to : 

2D''^ - D' coth(Ka) - \ sech^{K,a) 



B 



8t'^ sinh^(Ka) 



(11) 



(12) 



In Fig. [6l we plot the logarithm of the prefactor 
for both exact numerical results and for the approxima- 
tion given above, as a function of D/a. The agreement 
is good, given the fact that we neglect the deformation 
of the barrier due to the external bias potential in our 
approximation, and the qualitative trend is correct. Fig. 
[7] shows the prefactor as a function of Fermi level, and 
again we find qualitative agreement. The prefactor is ex- 
ponentially sensitive to system parameters, due to the 
fact that it is proportional to J(0); the dependence of B 




FIG. 6: Logarithm of quadratic bias-dependence prefactor 
versus spacer thickness D. Parameters are U = 1.75t, -Ef = 
—OAt, npi = np2 = 7. 
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FIG. 7: Logarithm of quadratic bias-dependence prefactor 
versus Fermi energy. Parameters are U — 1.75t,D/a = 
7, riFi = nF2 = 7. 



on system parameters is much weaker. For all parame- 
ters we've considered, we find that the zero bias exchange 
coupling between layers is antiferromagnetic in the tun- 
neling regime, and that the application of a bias increases 
this antiferromagnetic coupling (i.e., the sign of the V'^ 
prefactor is the same as the sign of J(0)). 

Recent experiments on MgO indicate that at bias po- 
tentials on the order of 0.5 V, the in-plane and out-of- 
plane torques are of similar magnitude, in stark contrast 
to metallic systems in which Fy ^ Fj,. The origin of 
this difference is the reduced area in A:-space for the in- 
plane torque in Fe-MgO. In metals all of the states in 
the Brouillon zone contribute to the in-plane torque and 
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only states near extremal points of the Fermi surface con- 
tribute to the out-of-plane torque, while in Fe-MgO only 
states near the F contribution to both in-plane and out- 
of-plane torques. 



IV. CONCLUSION 

In this work, we have derived a general expression for 
magnetic exchange interactions that is valid in nonequi- 
librium systems in the ballistic regime. We applied this 
formula to magnetic spin valves to study how the in- 
terlayer exchange coupling depends on an external bias 
voltage. We find that for metallic systems in the small 
confinement limit, the interlayer exchange coupling is an 
oscillatory function of the bias, but that the period of the 
oscillations is so long that seeing a change in interlayer ex- 
change coupling is not experimentally feasible. For tun- 
neling systems, and for any system in which electrons 
emanating from the L and R leads contribute equally to 
the interlayer exchange coupling and there is a symmetric 
bias potential drop, we find that the interlayer exchange 
coupling depends quadratically on the bias voltage, which 
is in agreement with recent experiments. 

We emphasize that the results we obtain are for sys- 
tems in the ballistic regime. The behavior of T±{V) for 
systems in the diffusive regime is generally different. In 
particular, a nonzero F^ at zero bias is a phase coherent 
effect; diffusive systems have a vanishing equilibrium Fj^. 
The oscillatory F-dependence described for metallic sys- 
tems is derived from modulating the wave-vector of co- 
herent states, and will not be seen for diffusive systems. 
Diffusive systems under bias generally have a nonzero 
F^. For symmetric systems, this dependence is quadratic 
and is therefore omitted from purely linear theories such 
as the spin circuit theory?*^. For asymmetric structures, 
there is a linear contribution tor±{V), which is described 
by linear, diffusive theories, where for example in the spin 
circuit theory it is related to the imaginary part of the 
mixing conductance 

There are other systems in which bias-altered exchange 
interactions can have important effects. It is possible, 
for example, that in single ferromagnetic layers under 
bias, intralayer exchange interactions can be affected by 
this bias, which would be manifest in a change of spin 
wave dispersion or Curie temperature. We are currently 
investigating these types of systems. 



APPENDIX: EXPANSION OF 
NONEQUILIBRIUM DENSITY MATRIX IN 
MAGNETIC ROTATIONS 

Here we describe in more detail the calculation of ex- 
change interactions in systems under bias. Our approach 
is based on the nonequilibrium Green's function method. 
In this method, we assume there is some central region, 
described by a Hamiltonian Hq, embedded between two 



semi-infinite reservoirs or leads. The left and right leads 
are described by the semi-infinite Hamiltonians Hj^ and 
i?R, respectively. We denote the coupling between the 
central region and the left /right leads by tl/tji. The 
retarded Green's function for the isolated central region 
is simply Gq{E) — [E + ii-j — Hq)^^ . The presence of 
the left lead is accounted for with self-energy Tji^{E) = 



is the so-called 



TLffEVt, where g£0 = {E + ir^ - H^^)-}, 
surface Green's function. This is the projection of the 
full semi-infinite lead Green's function onto the subset of 
lead sites which couple to the central region (for the one- 
dimensional model we consider, this subset is denoted as 
the (0,0) element of the full left lead Green's function). 
The self-energy for the right lead is defined similarly. 
The total retarded Green's function for the center region 
plus leads is then C'^E) = [E + irj - Hq - T^l - Ej;^)"\ 
In the following we combine left and right self energies: 

The so called "lesser" Green's function is given by the 

r<< — 1 /^rv</^o31i32 

The contribution to 



well-known Keldysh equation 

(- 



where S^^p 



^L/R 



the density matrix from states emanating from the L/R 
lead is: 



Gt/^{E) dE (A.l) 

-OO 

(A.2) 



Pl/r 



= ^ / G^^t/KG'-dE 



We consider how the density matrix changes upon the 
addition of some perturbation H' to the Hamiltonian. 
To this end we consider how the integrand of Eq. (jA.2p 
changes, and omit the energy integral for clarity. The 
integrand is given more explicitly as: 



1 



27r \E- He -S'' 



I 



E-Hc-Y'' 



(A.3) 



We restrict our attention to perturbations H' which are 
only present in the central region; H' leads to a new 
Green's function Gf , given by: 



Gf - ^ 



1 



1 



27r \ E - He- H' -Y."- 



1 



Hc-H' 



{kA) 



To first order in iJ', Gf is: 



Gf w — G"'(l + i7'G'')E<G"(l + i7'G") 

K G< + G'^H'G< + G<H'G''. (A.5) 

All of these objects are represented in a real-space 
basis, and we consider perturbations about a coUinear 
ground state. The components of the unperturbed 
Green's functions are then labeled with site index 
and are diagonal in spin space: 



Gi,j — 



G 





G^tj 



(A.6) 
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In all 2 X 2 matrices, we indicate explicitly site labels 
with indices and also spin label. As before, we define 
Ai.j = [H] — Hi)^j. We consider tilting some por- 
tion (sites {«'}) of the magnetization by a small angle 6, 
so that H' is given by: 



E 

ie{i'} 



A,: 



cos(6') - 1 sin(6') 
sin(6') 1 - cos(6') 



(A.7) 



There can in principle be spin-dependent terms that are 
off-diagonal in site index, but these are usually much 
smaller than site-diagonal terms and omitted here. In- 
serting Eq. (IA.7P into Eq. (|A.5P yields: 

G< -G< = 



GUcoa{9) 



1)A'G< 



sin(6l)A'G< 



Gl sm{0)A'Gf G'[(l - cos(e'))A'G< 



h.c. 



(A.J 

Where for notational clarity we've introduced matrices 
A', where AJ ^ = A^ ^ if i G {i'}, and otherwise. In 
the above h.c. stands for "Hermitian conjugate" . From 
the perturbed density matrix we find the ensuing spin 
densities (to first order in 9): 

ml. = 29 Re [(G^ A'G< + GTA'Gf) 



261 Ini 



(GfA'G 



GlA'Gf] 



2J 



(A.9) 



As discussed in the introduction, the exchange torques 
in a multilayer can be found by computing the out-of- 
plane torque. These torques result from the misalign- 
ment between the spin-dependent Hamiltonian and the 
spin density. To find this misalignment, we need the vec- 
tor components of the spin-dependent perturbed Hamil- 
tonian H + H', given as: A^ = Tr [{H + H') aa] , {a = 
x,y,z}, where aa are Pauli spin matrices. If the layers' 
magnetizations span the x — z plane, then the exchange 
torque is: 



r„ = A, 



A, 



(A.IO) 



Expanding to first order in 9 leads to the following expres- 
sion for the coupling due to states from the L/R leads: 



L/R 



Re 



Tr 



A'G| A'G]^L/R 



A'qA'G<L/R 



dE . 
(A.ll) 



In equilibrium, the total lesser Green's function G* 



- G^ can be written in terms of the retarded Green's 



function alone: G^ = (G'" — G"), in which case the 
expression for the total exchange torque becomes: 

1 f^^ 

J = -Im / Tr [A' (G^ - GJ) - A'Gf A'G^^] dE . 

^ J — OO 

(A.12) 



Combining Eq. (|A.12[) for energies in which both leads 
are occupied, and Eq. (|A.lip for energies in which only 
one lead is occupied yields Eq. 

We should note that in using the non-self-consistent 
spin density from a system rotated slightly out of 
coUinearity to find torques, we are implicitly making use 
of the magnetic force theorem. In addition, we utilize the 
same assumptions as Ref. S, namely that the orientation 
of spin magnetic moments on each lattice site is a well- 
defined quantity, and that these spin orientation degrees 
of freedom are "slow" with respect to the electronic de- 
grees of freedom (adiabatic approximation) . The validity 
of these approximations have been considered in detail in 
previous works^^. 

We make some final comments to compare our re- 
sult for exchange interactions with previous equilibrium 
expressions for exchange interactions. First, in similar 
spirit to Ref. Is!, we can attempt to define the pairwise in- 
teraction of two spins Ji J by finding the exchange torque 
present on atom i when i and j are both rotated, and 
subtracting off the torques present on i when i and j 
are rotated separately. Under nonequilibrium conditions 
however, the following sum rule is violated: 



Ji 7^ Jij ■ 



(A.13) 



implying that a pairwise interaction is not properly de- 
fined in nonequilibrium conditions. This is again due to 
the fact that there is a net flux of angular momentum 
into the system from the current-carrying electrons. 

Finally we note that the approach taken in Ref. d 
to finding exchange interactions relies on calculating the 
change in single particle energy AE upon rotating a por- 
tion of the magnetism away from coUinearity. It can be 
shown that in equilibrium, this change in energy can be 
expressed in the language of Greens functions as: 



AE 



1 



Tr [H' (p' + p)] 



(A. 14) 



where p' is the density matrix corresponding to the per- 
turbed lesser Green's function Gf. We have found that 
this expression works in nonequilibrium conditions as 
well, and gives identical results to the exchange torques 
found using spin densities. 



1. In-plane torques 

With the expansion of the nonequilibrium density ma- 
trix given by Eqn. (|A.5P , it is a simple matter to find the 
spin transfer torque, which we do here for completeness. 
As in Ref. il^, we find the out-of-plane spin density to 
determine the spin transfer torque present on a subsys- 
tem of atoms {j}, when atoms {i'} are rotated by a small 
angle 9. Assuming the two layers span the x — z plane, 
the y-component of the spin density is called for: the out- 
of-plane spin density in a spin valve determines the spin 
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IP 



\A\m„ 



torques. The spin transfer torque is thus F 
As before, we use the expansion of m given by Eq. (|A.9[) 
to obtain: 



^ip 



sm{0) 



Im 

h 



Tr 



A'G^A"G<- 



A'G"j;A"G< 



dE 



(A.15) 



Where the trace is over orbitals of the subsystem {j} 
in question, A' and A" are the spin-dependent Hamilto- 
nian for orbital sets {j} and {i}, respectively. Eq. (jA.lSp 
has the utility of expressing in plane torques in terms of 
coUinear objects, simplifying the numerical task of calcu- 
lating these torques from first principles relative to cal- 
culating full noncoUinear spin configurations, as in Refs. 
[H and|l7|. 



of a single potential Va,b is described by the reflection 
and transmission coefficients ta^b and Ia^b- We suppose 
the potential to be confined to sites i within a range 
i € {i'}, and that the "left-most" site is i^. For our 
purposes, it suffices to consider the change in Green's 
function SGij induced by the potential for indices i,j on 
the same "side" of the potential, i.e. i,j < il. Then, for 
small Va (see Fig. dH])): 

SGi,j = {GqTaGq)^- 

= ihvkGo^^ ^^/A Goti^.^) . (A. 21) 



OOOOdODOOOOO 



2. Bias dependence of interlayer exchange coupling 

Here we take the limit of large spacer thickness and 
small reflectivities, and reformulate J in the language of 
Refs. and0. We consider first the equilibrium case as a 
precursor to the nonequilibrium case. We can rewrite the 
Green's functions in terms of free space Green's functions 
Go: 

G'^ = GI + GI{Ta + Tb)GI , (A. 16) 



Ta,b = 



g«fc|i-i' 
ihvk ' 

1 - GIVa,b 



(A.17) 
(A.18) 



T includes the influence of the potentials Ya,b exactly. 
When the potentials V is small, then T V and G is 
approximately given as: 



G-- = GI + GI{Va^Vb)GI 



(A.19) 



Inserting Eq. (|A.19p into Eq. (jA.12|) and retaining the 
lowest order terms in V leads to a useful form of J: 



J 



1. 



-Im Tr 



yi-yilGo {v\~v)^\G^dE 



(A.20) 



In the spirit of previous works on interlayer exchange 
coupling, we next rewrite Eq. (|A.20p in terms of spin- 
dependent reflectivities of the magnetic layers. The effect 



FIG. 8: Schematic for effect of potential V on Gij: the path 
from site i to site j can contain a reflection r off of the left 
edge of the potential V. 



Using the fact that Va,b ~ applying Eq. (IA.21|) 

twice, and the form for Gq given in Eq. (jA.17p leads to: 



J = -Im 

TT 



(A.22) 



where D — a(ns + 1). For the nonequilibrium case, we 
insert Eq. (IA.I9j) into Eq. (jA.lip . and keep lowest order 
in V: 

Jl,r = hm I""" (yj - V^) {gITlmG^, (vl - Vi) Gg 
+GS (vl - vl) GSFi.flGg) dE . (A.23) 



In terms of reflectivities and phase factors, this becomes: 



Jl = , 
Jr = — Im 

TT 



B ' B 



2ikD 



dE . 
(A.24) 



Evaluating the integral in the large D limit leads to Eq. 
dSD of the text. 
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